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Overview 

This report is the official research summary concerning Cooperative Agreement 
NCC- 1-382 between the Pennsylvania State University and NASA Langley Research 
Center. Although the cooperation between Penn State and NASA Langley is continuing 
beyond the official dates of performance, this summary report documents the status of 
research as of March 2002 and gives a brief, comprehensive summary of significant 
accomplishments during the duration of the cooperative agreement. 

Objectives 

The main proposed tasks of Cooperative Agreement NCC- 1-382 were: 

1. developing MSC/N ASTRAN DMAP language scripts to implement the Soize fuzzy 
structures approach for modeling the dynamics of complex structures, 

2. benchmarking the results of the new code to those for a cantilevered beam in the 
literature, and 

3. testing and validating the new code by comparing the fuzzy structures results to NASA 
Langley experimental and conventional finite element results for two model test 
structures representative of aircraft fuselage sidewall construction: 

3 A. a small aluminum test panel (SLP, single longeron panel) with a single longitudinal 
stringer attached with bolts, and 

3B. a 47 by 72 inch flat aluminum fuselage panel (AFP, aluminum fuselage panel) 
including six longitudinal stringers and four frame stiffeners attached with rivets. 

Accomplishments 

As of March 2002, all the major objectives of the cooperative agreement were 
met, although additional work needs to be performed for objectives 3A and 3B above and 
additional documentation needs to be produced to fully transfer the fuzzy structures 
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techniques for in-house use at NASA Langley. The accomplishments under each 
objective will now be reviewed. 

Regarding Objective 1, the MSC/NASTRAN DMAP codes were written during 
the spring and early to mid summer of 2000. Approximately 800 lines of DMAP (direct 
matrix abstraction program) statements are involved. The DMAP codes are accompanied 
by a MATLAB script to generate parameters needed by the fuzzy structures analysis 
predictions. 

Regarding Objective 2, during the late summer and fall of 2000 the new DMAP 
codes were applied to a test problem published by C. Soize (a cantilevered beam) for 
which fuzzy structures results were available. The new DMAP results and those 
previously published by Soize agreed very well, and this work was presented at a spring 
2001 AIAA meeting, see Attachment 1. Hence, the DMAP code benchmarked very well, 
and the research then proceeded with Objectives 3A and 3B. 

The single longeron panel (SLP) modeling took place during the spring and early 
summer of 2001. The SLP is described in detail in Attachment 1. These results were 
presented at NASA Langley and at the spring 2001 meeting of the Acoustical Society of 
America. Here the nut and bolt joints were assumed to have uncertainty. For a 
representative comparison plot of the drive point frequency response, see Attachment 2. 
In general, the mean fuzzy structures analysis result is overly damped compared to the 
measured response. It was noted that this structure had very little complexity and very 
little internal damping. Hence, in retrospect, this probably was not the best candidate 
structure on which to apply the fuzzy structures analysis procedures. 

The aluminum fuselage panel (AFP) fuzzy structures analyses were performed 
during the late summer and fall of 2001 and continue into the spring of 2002. The 
fuselage panel is described in detail in Attachment 1. For this panel uncertainty is 
assumed in the rivet attachments of the frame stiffeners. Preliminary comparisons 
between measured data, conventional finite element results, and fuzzy structures finite 
element results are encouraging. See Attachment 3, again showing a drive point 
frequency response. Obviously for this large model (with 60,000 degrees of freedom) the 
fuzzy structures analysis procedures behaved pretty much as predicted. It is now clear 
that the fuzzy structures DMAPs will produce the type of results that will be useful in the 
uncertainty modeling of large dynamic structures of interest for studying aircraft interior 
noise. 

Additional work regarding the AFP is ongoing, and this research will be written 
up and submitted for publication in a peer-reviewed journal. 


Attachments 

1. V. W. Sparrow and R. D. Buehrle, “Fuzzy structures analysis of aircraft panels in 
NASTRAN,” AIAA-2001-1320 presented at the 2001 AIAA Structures, Structural 
Dynamics, and Materials Conference, Non-Deterministic Approaches Forum, Seattle, 
WA, April 16, 2001. 

2. Representative frequency response comparison plot showing the application of the 
fuzzy structures analysis procedure to the single longeron panel (SLP). Plot taken 
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from the presentation of V. W. Sparrow and R. D. Buehrle, at the 141 st Acoustical 
Society of America Meeting, Chicago, IL June 6, 2001, “Application of fuzzy 
structures analysis to aircraft panels,” J. Acoust. Soc. Am. 109 (5, Pt. 2) 2410 (2001). 

3. Representative frequency response comparison plot showing the application of the 
fuzzy structures analysis procedure to the aluminum fuselage panel (AFP). 
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ATTACHMENT 1 


AIAA— 2001— 1320 


FUZZY STRUCTURES ANALYSIS OF AIRCRAFT PANELS IN NASTRAN 

Victor W. Sparrow* 

The Pennsylvania State University, University Park, PA 16802 
Ralph D. Buehrle + 

NASA Langley Research Center, Hampton, VA 23681 


Abstra ct 

This paper concerns an application of the fuzzy 
structures analysis (FSA) procedures of Soize to 
prototypical aerospace panels in MSC/NASTRAN, 
a large commercial finite element program. A 
brief introduction to the FSA procedures is first 
provided. The implementation of the FSA methods 
is then disclosed, and the method is validated 
by comparison to published results for the forced 
vibrations of a fuzzy beam. The results of the 
new implementation show excellent agreement to 
the benchmark results. The ongoing effort at 
NASA Langley and Penn State to apply these 
fuzzy structures analysis procedures to real aircraft 
panels is then described. 

Introduction 

In 1986 Christian Soize was at the Office Na- 
tional d 'Etudes et de Reoherches Aerospatiales 
(ONERA) 1 * in Chatillon, France and introduced the 
concept of fuzzy structures 2,3,4,5,6 in the dynamic 
analysis of large scale structures. Having nothing 
to do with the use of fuzzy sets or fuzzy logic, 
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a fuzzy structure is composed of (1) a well char- 
acterized master structure and (2) attached fuzzy 
substructures whose details are either unknown or 
are imprecisely known. In the context of a tradi- 
tional finite element model, the master structure 
would be that part of the structure amenable to 
conventional FEs. The fuzzy substructure would 
be those parts of the structure inaccessible to tra- 
ditional modeling: missing degrees of freedom, fine 
scale structural details, attached equipment, etc. 
An implementation of the fuzzy structures concepts 
in FEs as originally conceived by Soize is called 
Fuzzy Structures Analysis (FSA), a term coined by 
Ruckman. 7 

Although there have been some attempts to 
implement alternative fuzzy structure approaches 
in the context of analytical modeling, 8,9,10,11,12 the 
original FSA methodology of Soize has great attrac- 
tiveness since it can be performed using existing FE 
models: no new special elements need be generated. 
What has been lacking, however, is an implementa- 
tion of FSA in a mainstream FE code. 

This paper describes an original implementa- 
tion of Soize’s FSA in the well known commercial 
code MSC/NASTRAN, 13 specifically its solution 
sequence 108, the dynamic frequency response mod- 
ule. Using NASTRAN’s Direct Matrix Abstraction 
Program (DMAP) environment and a MATLAB 14 
module to provide appropriate input parameters, 
the new implementation allows one to include fuzzy 
substructures along with any structural model al- 
ready existing in NASTRAN. The new DMAPs 
have been benchmarked against published results 
for vibrations of beams with attached fuzzy sub- 
structures, and these benchmarks are reported here 
for the first time. This paper also describes an on- 
going application of the new FSA implementation 
to two representative aircraft fuselage panels under 
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study at the NASA Langley Research Center. 

This paper is arranged with the following sec- 
tions: Firstly, the background of Soize’s fuzzy 

structures analysis is provided, focusing on both 
the overall concepts of the method and the required 
mathematical formulation. Secondly, this paper 
describes the specific implementation of the fuzzy 
structures analysis procedures in MSC/NASTRAN 
and MATLAB. Next a benchmark solution is de- 
scribed that was used to test the procedures. Two 
aircraft panels are then described which are the 
focus of our ongoing study in applying the fuzzy 
structures analysis procedures to real aircraft pan- 
els. 

Note that the acronym FSA is used for fuzzy 
structures analysis throughout this paper. Also 
for consistency with Soize’s works, the e +twt time 
convention is used. 

Background of FSA 

Motivation 

Soize introduced fuzzy structures for modeling 
systems with a huge number of degrees of freedom 
with a reduced size model. A good example of a 
fuzzy structure, a system too complicated for con- 
ventional modeling, would be the electronics under 
the cockpit of a commercial airliner. This mechan- 
ical system cannot be modeled directly because it 
is so complex. However, it must be modeled some- 
how since it contributes to the overall dynamics and 
structural acoustic response of the airplane. This 
particular fuzzy structure is important for interior 
cabin noise upon landing since it may be directly 
above or near the forward landing gear of the air- 
liner, and the landing gear vibrations can travel 
along the fuselage and cause acoustic radiation into 
the cabin. 

Another aim for FSA procedures is to account 
for uncertainty of small scales in structural dy- 
namics models. One source of uncertainty is in the 
variability in the construction of large airliners. An- 
other source of uncertainty is in the finite element 
modeling procedure itself, since no individual can 
possibly explicitly incorporate into a deterministic 
finite element model all of the fine structure in a re- 
alistic airliner model. One idea put forth by Jeffrey 
Cipolla 15 in 1995 was to think of FSA as a large 
eddy simulation (LES) for large dynamics prob- 
lems in analogy with how fluid dynamicists handle 
fluid scales too small for explicit representation in 


a direct numerical simulation (DNS). 

The computational cost is low for FSA com- 
pared with the conventional Monte Carlo modeling 
of uncertainties in finite element simulations. In 
1986 Soize showed one could achieve similar results 
between a FSA and a traditional Monte Carlo model 
of an uncertain structure but with many fewer FE 
solves. 2 * This computational savings, along with 
keeping only a minimum number of DOFs in the 
FE model, is the great advantage of FSA. 

Previous work 

The algorithm for implementing Soize’s FSA 
for complicated structures has been available in the 
literature for many years. A somewhat complicated 
algorithm to be described briefly below, the proce- 
dure involves making a few runs of a conventional 
structural dynamics finite element analysis program 
with a large number of carefully constructed right 
hand sides (forcing functions). The results obtained 
using each of the forcing functions is then handled 
using standard statistical techniques. 

In 1994 the Carderock Division of the Naval 
Surface Warfare Center tasked one engineer, Chris 
Ruckman, to travel to France to investigate Soize’s 
fuzzy structures methods in detail. One of Ruck- 
man’s goals was to determine the feasibility of 
implementing the fuzzy structures procedures in 
NASTRAN, and he wrote a short report 16 on how 
he thought one could accomplish this. He also 
published an even shortened version of the report 
as a paper at conference. 7 After the work was 
written up, changing priorities and reduced bud- 
gets at NSWC and the Office of Naval Research 
stopped further development of implementing the 
fuzzy structures procedures in NASTRAN. Ruck- 
man left NSWC, and he moved on to other areas. 
However, Ruckman correctly predicted that a fuzzy 
structures analysis in NASTRAN was possible. 

At about the same time the company now 
known as BBN Technologies of Cambridge, MA 
also was tasked by the Office of Naval Research 
to include a fuzzy finite element capability in their 
program SARA, a coupled fluid and structure finite 
element solver. 17 This development was accom- 
plished in 1994 using the Pierce/Sparrow/Russell 8 
and Strasberg/Feit 12 fuzzy structure models instead 
of Soize’s model. Because of reduced budgets, how- 
ever, no detailed testing or validation of the BBN 
implementation, to the best of our knowledge, was 
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ever made. SARA is not used in the commer- 
cial aircraft industry. However, since NASTRAN 
is more widely in use, the present authors decided 
that the prudent course was to implement Soize’s 
procedures in NASTRAN, putting them into the 
mainstream for testing and validation by NASA 
Langley, industry, and university researchers. 

Philosophy 

The general philosophy of FSA is to model 
as much of structure/fluid with conventional finite 
elements as possible. This portion of the structure, 
accessible to conventional modeling techniques is 
called the master structure. To add the influence 
of additional fine detail or uncertain attachments 
to the master structure, one modifies the boundary 
conditions of the master structure. See Fig. 1, where 
the Z's indicate boundary impedances affecting 
the master structure. This addition is called the 
fuzzy substructure. By only incorporating the 
details of the fine or uncertain structure through 
the boundary conditions only, the total number 
of degrees of freedom (DOFs) in the model is 
unchanged. This helps keep the size of the finite 
element model from increasing to a point where 
conventional FE solvers are unusable. 

It should be noted that the fuzzy substructure, 
in the present analysis, only affects the master 
structure through the FE nodes. One only places 
the fuzzy substructure at nodes where one expects 
the uncertainty or fine structure to exist. 

Fuzzy impedances 

Each additional boundary condition in the 
fuzzy substructure is modeled as the impedance 
of a one degree of freedom attachment. See Fig. 2 
where the mass is connected to the master struc- 
ture through a stiffness K and a damping constant 
C and where u) p = (K//i)s is the natural frequency. 

Here £ = 2 ^ is the critical damping constant as- 
sumed to be small, (0 < f < 1). Assuming the 
support has a normal velocity x = v = ve wt with 
driving frequency cj, one can easily show 

p -fjul - 1 + 4£ 2 ) -I- iuj (2^£Wp) 

2 = _ — - A_ L (1) 

c +- (£)«’] 

for the impedance of the attachment. 

To make the oscillator fuzzy /uncertain the 
following assumptions are made about the mass, 


damping, and modal density n of the oscillator: 

/i — /x (1 -f Fi) 

£ = £(l+*2) (2) 

n = n (1 -b Y3) 

where _ indicates mean value, and where the Y { 
axe independent random variables, in practice small 
fractions of unity. 18 

A one DOF oscillator having modal density 
may be unfamiliar, but it simply is equivalent to the 
oscillator having uncertainty in natural frequency 

OJp. 

Incorporation into the dynamic model 

In the following paragraphs U indicates a col- 
umn vector of generalized displacements, F a col- 
umn vector of generalized forces, and square brack- 
ets [ ] indicate matrices. A conventional structural 
dynamics vibration problem would be formulated 
for a particular angular frequency u) using the FE 
equation 

iw ([-^master structure]) U = F (w) . (3) 

Adding the effect of the fuzzy substructure consists 
of modifying the dynamical equation at appropriate 
nodes as 

iiV ([^master structure ] + [Zfu»y])tf = . (4) 

Decomposing the fuzzy substructure into its mean 
and random components one has [Zf uzzy ] = [Z] + 
[^rand] yielding 

IUJ ([Zdet] + [Zrand]) U = F (u>) (5) 

for the dynamical system, where [Zdet] — 

[^master structure] ”1" l^j* Soize then pulls Out [-Zdetjj 

iu[Z det \([I}-[T])U = F(u) , (6) 

an equation which can be solved by a perturbation 
analysis. Here [T] = — [^det] _1 [^rand] has entries 
assumed small compared to unity. Soize defines 

M 

[Z rand ]=]T X< [Zi] (7) 

i= 1 

where the [Z t ] are the contributions to [Z ra n d] 
by each random aspect of the fuzzy. There is 
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only 1 non-zero element in each [Z x \. There are 
M = L x J x d contributions to the sum: 

L = number of attachments 
J — number of directions (3 in general) 
d = number of fuzzy parameters (usually 3 4) 
and the \i are independent, centered, uniform 
random variables. Note that the construction of 
the [Zi], too lengthy to give here, is completely 
described in the works of Soize 2,3,4,6 as are other 
details concerning the solution procedure. 

Formally, if it could be constructed, the solu- 
tion of Eq. (6) is 

U = {{I}-\T])- l (iu\Z <iet \r l F( ! j) . (8) 

Because ||T|| < 1 one can write 


3. Postprocessing of results gives mean and vari- 
ance of response. 

Here are additional details: 

Step 0: k = 0. The initial step gives by solving 
iw[Z det ] Z7 (0) =F(w) . (13) 

This is equivalent to solving Eq. (10) omitting the 
sum completely. Note that the solution 
of Eq. (13) is completely deterministic. 

Step 1: k = 1. Get for 1 < i < M by solving 
iuj[Z det ] Z7' 1} = -iaj[Z t }U {0) . (14) 


([I]-[T])- 1 = [I} + J2^ k ■ ( 9 ) 

fc= 1 

So Eq. (8) becomes 

(io) 

*=i / 

where 

M 

[T] = -[Z det )- 1 Y,Xi{Zi] ■ (11) 

v — 1 

The solution is constructed by expanding out 
Eq. (10): 

M MM 

V = u m + £ +ZE + • • • 

i=l n = L 1 2 — 1 

( 12 ) 

Each of the U are solved for easily, and each 
of these solves is completely deterministic. The 
are the only random variables in Eq. (12). 

Solution sequence 

The usual solution sequence follows from solv- 
ing Eq. (10) with increasing; orders of the index k 
in these steps: 

0. FE solve with master structure and mean fuzzy 
only. 

1. Another FE solve of master structure with 
specially constructed right hand sides (RHSs) 
using results of 1st solve (Step 0). 

2. If necessary, another FE) solve of master struc- 
ture using specially constructed RHSs using 
results of 1st and 2nd solves (Steps 0 and 1). 



This is equivalent to solving Eq. (10) keeping only 
the first term of the surn - Usually 

Eq. (14) is implemented using 1 solve with M 
RHSs for each frequency. Each frequency must 
be solved separately since [Z^et] turns out to be 
frequency dependent. Further, this k = 1 solve, 
again completely deterministic, usually is enough 
to represent Eq. (12) with sufficient accuracy. 

Step 2: k = 2. If necessary one further solves for 
U^\ 2 for 1 < i\ < M and 1 < 12 < M by solving 

iu\Z^}U™ 2 = -i“[Zi 2 )U { " . (15) 

This is equivalent to solving Eq. (10) keeping the 
first two terms of the sum. Here Eq. (15) 

is implemented using 1 deterministic solve with 
M 2 RHSs. Since a solve with M 2 RHSs could be 
prohibitively expensive, one hopes that this step is 
not necessary. 

Step 3. Postprocessing: 

Now that the u[ l \ and if necessary 
components of Eq. (12) have been found, one can 
find the statistics of the solution. Ruckman 7,16 
elucidated this procedure. 

Clearly the solution is in the form of gener- 
alized displacements and the expectation (average) 

of U is simply obtained in step 0 of the so- 
lution sequence. However, the expectation of any 
magnitude squared value is 


E{\r\ 2 } 



i= 1 


(16) 
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for k = 1 when r is any scalar value such as 
displacement, velocity, acceleration, pressure, etc. 
The variance is further given for k — 1 as 

M 9 ^ M , 4 

-|£ r!"’ 

1=1 1=1 

MM 9 

+ 2^S[Re{r«V>}] . 

i=l j — 1 

(17) 

Implementation 

Overview 

To implement the solution procedure just de- 
scribed it was decided to use two consecutive NAS- 
TRAN DMAP runs, including modifications to the 
SOLUTION 108 (direct frequency response) mod- 
ule. The NASTRAN database capabilities would be 
used to minimize matrix recalculations between the 
two runs. The first NASTRAN DMAP run saves all 
the necessary values through a simple execution of 
the usual 108 sequence, performs precalculations, 
and saves needed information to the NASTRAN 
database. The second NASTRAN DMAP run ac- 
tually executes the k = 0 and k — 1 steps in the 
solution sequence described above as well as per- 
forming the solution postprocessing. An extensive 
MATLAB script is executed after the first NAS- 
TRAN run to precalculate information for the [Zi] 
matrices, required for the second NASTRAN run. 

The first DMAP code consists of 123 lines of 
DMAP inserted into usual SOLUTION 108 (direct 
frequency response) bulk data file. Here is where 
the usual 108 sequence is used. It also saves the 
database and counts the number of fuzzy DOFs. 
The second DMAP code contains over 650 lines of 
DMAP. It uses the NASTRAN database saved in 
the first DMAP, but is essentially a stand alone 
analysis. 

The MATLAB code, “makefuzzy,” provides 
inputs for the second DMAP code since NASTRAN 
can’t perform 2-D numerical integration, required 
for generating the [ Z x ] matrices. 

Details on FSA DMAP 2 

Usually NASTRAN represents the mass, 
damping, and stiffness matrices straightforwardly 
as 

iui [Z det ] = -ui 2 [M] + iu> [.D] + [K\ . (18) 


However, in FSA the damping and mass matrices 
become frequency dependent as 

iw [Zdet] = ~ UJ 2 ([M] + [fifuz (w)]) 

+ iuj ([£> (w)] + [/fuz (w)l) + [ K (w)] 

(19) 

as they now include the real and imaginary con- 
tributions of the fuzzy substructure, [Rf uz {u)] and 
[/fuzM], respectively. 

Hence the FSA DMAP 2 must solve for one fre- 
quency at a time. Note that this is not NASTRAN’s 
default procedure. Usually MSC/NASTRAN at- 
tempts to perform the direct frequency response 
for all frequencies in solution 108 during the same 
solver pass. 

To explicitly clarify what has to occur, a 
solve for the 1st frequency would take the form 
zu;[Zde t(^i)] = F(uq), then one solves 

UJ [Z del (wi)] u[ 1] = -iw[Zi(u> i)]C7 (0) (u/i) (20) 

for each i. Similarly, the corresponding solve for the 

2nd frequency would be iu[Z ^{^ 2 )] = F(uj 2 ), 

then one solves 

iuj [Zdet (wa)] u\ l) = -iu) [Zi (wa)) U {0) (u> 2 ) 

again for each i, etc. 

User procedure 

When running a typical FSA, the following 
procedure is followed: 

• Run FSA DMAP 1 

- (use scr=no command line option) 

- Input: nodes designated as fuzzy in a SET 

- Output: .DBALL database file 

- Runtime: about same as SOLUTION 108 

• Run Matlab Code “makefuzzy” 

- Input: frequencies for analysis, 

fuzzy parameters 

- Output: file fuzzyforphasela2.txt, 

matrix element values 

- Runtime: 5 min on a DEC Alpha work- 
station for 100 frequencies. 

• Run FSA DMAP 2 

- Input: fuzzyforphasela2.txt, 

“FSA DMAP 1” .DBALL file 

- Output: mean squared displacement, etc 
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- Runtime: 2 x (SOLUTION 108 runtime) 
4- time for RHS calculations. 

The number of RHSs = 9 x number of fuzzy 
nodes. This is because in the present implementa- 
tion only displacement DOFs (3) are assumed to be 
fuzzy, and only 3 fuzzy variables are assumed for 
each DOF. The 3 fuzzy variables are mass, critical 
damping, and modal density given in Eq. (2). The 
amount of deviation for these variables from their 
mean values are specified by the user as dispersion 
parameters Ai,A 2 ,A 3 in the MATLAB procedure 
“makefuzzy. “ 

Testing and validation 

The resulting DMAP/Matlab procedures have 
been successfully applied to a benchmark prob- 
lem published by Soize, a fuzzy simply supported 
beam. 5 

Model parameters 

The base beam (master structure) parameters 
were taken directly from Soize’s work 5 . The beam 
mass was 20,000 kg, and the attached oscillators 
have a 0.002 damping rate and a mass of 1.8 kg. 
The dispersion parameters were matched to be 

X x = A 2 = A 3 = 0.001. 

However, the beam modeled here did have 
some differences to Soize’s beam. Soize’s model had 
2500 FE nodes and was modeled using a 0.05 Hz 
frequency resolution. He used 1213 1-DOF fuzzy 
attachments in the y (“up”) direction, with n = 
0.035, 0.043 for mean modal density on two distinct 
patches of oscillators. In the present work only a 
200 FE node model was employed using a 1 Hz 
frequency resolution from 1 to 100 Hz. Here 300 
1-DOF fuzzy attachments were distributed over the 
two patches specified in Soize’s problem, but they 
were in all 3 directions (re, y, and z) and they 
all used the same value, 0.035, for mean modal 
density. The 300 fuzzy attachments were located 
at 100 fuzzy nodes, corresponding to 900 RHSs in 
the k — 1 FE solution. In summary for the present 
FSA DMAP 2, 100 k = 1 individual FE solves were 
made, 1 solve for each frequency, with each solve 
containing the 900 RHSs. 

The run times on a NASA Langley Structural 
Acoustics Branch DEC Alpha workstation were 

• 1 minute — FSA DMAP 1 

• 5 minutes — “makefuzzy “ Matlab file 

• 6 h, 7 minutes — FSA DMAP 2 


Clearly the codes are not yet optimized for 
speed, but these numbers provide some evidence of 
the substantial time required for the procedure. 

Results 

The plot shown in Fig. 3 clearly indicates the 
similarity of the results for the mean squared ac- 
celeration at a point on the beam published as Fig. 
7(a) of Soize. 5 Here the solid line gives the fuzzy 
structures result and the dashed line gives the re- 
sult for the master structure alone. (The additional 
thin solid line in Soize’s Fig. 7(a) is one realiza- 
tion of a Monte Carlo simulation for comparison.) 
Sparrow’s beam model has a coarser discretization 
(200 FE nodes) and frequency resolution compared 
to Soize’s model (2500 FE nodes), but otherwise 
the plots agree very well in magnitude. A close in- 
spection of the two results shows that there is a 
slight discrepancy in the location of the resonances, 
and this again is due to the different discretizations 
of the models. These results clearly indicate that 
these new fuzzy structures procedures in DMAP 
are working correctly and are ready for application 
to the more realistic panel structures of interest to 
the aerospace and transportation industries. 

Note also in Fig. 3 the greatly constricted 
fluctuations of the structure’s acceleration with the 
addition of the fuzzy substructure. In essence, 
the fuzzy substructure acts to damp the master 
structure vibrations. 

Prototypical aircraft panels 

The present research at Penn State and NASA 
Langley is now focused on applying the new FSA 
procedures to realistic, prototypical aircraft panels. 
The investigation is centered on modeling the junc- 
tions between the longerons and skin and frames 
and skin in such panels as fuzzy /uncertain due to 
the inherent impreciseness of such joints. Note that 
this is a departure from Soize’s original applica- 
tion of the FSA procedures for uncertain attached 
equipment in large dynamical systems. 

Single longeron panel 

The first panel consists of a .04 inch thick 
aluminum skin with a single longitudinal stringer 
dividing the panel skin into two equal bays of 37.75 
by 7.625 inches. A single line of thirty 6-32 x 
3/8 screws attach the longeron to the panel skin. 
The panel is mounted in a rigid frame to simulate 
clamped boundary conditions. It is referred to as 
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either the single longeron panel, or the two-bay 
panel. 

Flat fuselage panel 

This second panel is representative of current 
aircraft construction but was manufactured without 
curvature to simplify the experimental and analyti- 
cal modeling. It is constructed of a 47- by 72-inch 
aluminum panel with a 0.050-inch skin with six 
equally spaced longerons and four equally spaced 
frame stiffeners. Single lines of rivets attach the 
stringers and frames to the skin. A bay is defined 
as a section of the panel skin that is bounded by 
the stringers and frames. The bay responses are 
the focus of the fuselage panel correlation efforts. 
For both the experimental and analytical studies 
on this panel, free-free boundary conditions were 
simulated. 

Progress in modeling 

Conventional FE modeling of the fuselage 
panel 19,20 has focused on different stiffener to skin 
attachment models and their effects on the pre- 
dicted bay motions. To characterize the panel 
dynamic response up to 1000 Hz, several methods 
of modeling the panel were examined. First, the 
required finite element mesh density of the panel 
skin was evaluated by performing a normal mode 
analysis of a single bay with clamped boundary 
conditions. The skin was modeled with linear plate 
elements. A mesh of 30 by 16 elements was found 
to provide a one-percent convergence on frequency 
and adequate spatial resolut ion to define the mode 
shapes through 1000 Hz. This resulted in 11682 
linear plate elements for the overall panel skin. 

For the FSA study, the stiffeners were modeled 
using one-dimensional beam elements with the ef- 
fective cross-sectional properties (area, inertia, and 
torsional constant) of the stiffeners. Offsets from 
the skin to the stiffener neutral axis are also in- 
cluded. The beam elements are created along a 
line consistent with the rivet line that attaches the 
skin to the stiffeners. The beam stiffener model 
of the panel contains approximately 60,000 degrees 
of freedom (DOFs). Similarly, the two-bay panel 
was modeled using a combination of linear plate 
elements for the skin and beam elements for the 
longeron. The resulting model of the two- bay panel 
has approximately 7900 DOFs. 

The FSA modeling effort is currently underway. 
Preliminary results indicate that one must be very 


careful in picking appropriate values of fuzzy mass, 
fuzzy damping, and fuzzy modal density as input 
parameters for the FSA procedures. If these studies 
are successful, plans are underway to apply the 
FSA procedures to more complicated prototypical 
structures including NASA Langley’s Aluminum 
Testbed Cylinder (ATC). 21 

Conclusions 

This paper has described an implementation of 
Soize’s fuzzy structures analysis (FSA) procedures 
in MSC/NASTRAN and MATLAB. The proce- 
dures have been validated by comparison to the 
results for a simply supported fuzzy beam already 
in the literature. Previous studies have shown that 
with FSA one can obtain similar results to a Monte 
Carlo simulation of uncertain attachments with a 
greatly reduced computational expense. The work 
now underway is aimed at the fuzzy modeling of 
the skin attachments in two prototypical aircraft 
panels, and detailed results are expected soon. 

Acknowledgment 

Work supported by NASA Research Coopera- 
tive Agreement NCC-1-382. 

Disclaimer 

Any opinions, findings, conclusions, or recom- 
mendations expressed in this publication are those 
of the authors and do not necessarily reflect the 
views of the National Aeronautics and Space Ad- 
ministration. 

References 

1. See www.onera.fr. 

2. C. Soize, “Probabilistic structural modeling 
in linear dynamic analysis of complex me- 
chanical systems. I. Theoretical elements,” 
La recherche Aerospatiale No. 1986-5, 23-48 
(1986). 

3. F. Chabas, A. Desanti, and C. Soize, “Prob- 
abilistic structural modeling in linear dy- 
namic analysis of complex mechanical systems. 
II. Numerical analysis and applications,” La 
recherche Aerospatiale No. 1986-5, 49-67 
(1986). 

4. C. Soize, “A model and numerical method in 


7 


American Institute of Aeronautics and Astronautics 



the medium frequency range for vibroacous- 
tic predictions using the theory of structural 
fuzzy/ 5 J. Acoust. Soc. Am. 94, No. 2, Pt. 1, 
849-865 (1993). 

5. C. Soize, “Vibration damping in low-frequency 
range due to structural complexity: A model 
based on the theory of fuzzy structures and 
model parameters estimation,” Computers 
and Structures 58, No. 5, 901-915 (1995). 

6. R. Ohayon and C. Soize, Structural Acoustics 
and Vibration (Academic Press, San Diego, 
1998), particularly Chap. 15. 

7. C. E. Ruckman and D. Feit, “A tutorial 
on Soize’s method for stochastic modeling in 
structural acoustics (fuzzy structures analy- 
sis),” in DE-Vol. 84-2, Proceedings of the 1995 
ASME Design Engineering Technical Confer- 
ences, Vol. 3B: Acoustics, Vibrations and 
Rotating Machines, pp. 241 246. Edited by 
R. Daddazio, M. Ettouney, N. Abboud, et al. 
ASME, New York (1995). 

8. A. D. Pierce, D. Russell, and V. W. Sparrow, 
“Fundamental structural acoustic idealizations 
for structures with fuzzy internals,” ASME J. 
Vib. Acoust., 117 339-348 (1995). 

9. D. A. Russell and V. W. Sparrow, “Backscat- 
tering from a baffled finite plate strip with 
fuzzy attachments,” J. Acoust. Soc. Am., 98 
(3) 1527-1533 (1995). 

10. A. D. Pierce, “Resonant frequency distribu- 
tion of internal mass inferred from mechanical 
impedance matrices, with application to fuzzy 
structure theory,” ASME J. Vib. Acoust., 119 
324-333 (1997). 

11. R. L. Weaver, “Mean and mean-square re- 
sponses of a prototypical master/ fuzzy struc- 
ture,” J. Acoust. Soc. Am. 101, No. 3, 
1441-1449 (1997). 

12. M. Strasberg and D. Feit, “Vibration damping 
of large structures induced by attached small 
resonant structures,” J Acoust. Soc. Am. 99, 
No. 1, 335-344 (1996). 


13. MSC/NASTRAN, (MSC Software Corp., Los 
Angeles, CA.) See www.mscsoftware.com . 

14. MATLAB, Ver. 5.0. (Mathworks Inc., Natick, 
MA). See www.mathworks.com . 

15. J. Cipolla, private communication, 1995; for- 
merly with Naval Undersea Warfare Center, 
Division Newport; now with Hibbitt, Karls- 
son, & Sorensen, Inc., Newport, R.I.. 

16. C. E. Ruckman “Preliminary report on im- 
plementing fuzzy structures analysis in NAS- 
TRAN,” Carderock Division Naval Surface 
Warfare Center Rept. CDNSWC-SIG-94/142 
(October 1994). 

17. SARA (BBN Technologies, Cambridge, MA). 
See www.gte.com/ AboutGTE/ gto/ bbnt/ 
apts/ offerings/ sara.html . 

18. V. W. Sparrow, D. A. Russell, J. D. Rochat, 
“Implementation of discrete fuzzy structure 
models in Mathematical' Int. J. Num. Meth. 
Engr. 37, 3005-3014 (1994). 

19. G. A. Fleming, R. D. Buehrle, and O. L. 
Storaasli, “Modal Analysis of an Aircraft Fuse- 
lage Panel Using Experimental and Finite El- 
ement Techniques, 55 Proceedings of the 3rd 
International Conference on Vibration Mea- 
surements by Laser Techniques, Ancona, Italy, 
June 1998. 

20. R. D. Buehrle, G. A. Fleming, R. S. Pappa, 
and F. W. Grosveld, “Finite Element Model 
Development and Validation for Aircraft Fuse- 
lage Structures, 55 Proceedings of the 18th In- 
ternational Modal Analysis Conference, San 
Antonio, Texas, February 2000. 

21. F. W. Grosveld, “Structural normal mode 
analysis of the aluminum testbed cylinder 
(ATC),” AIAA Paper 98-1949, 39th Struc- 
tures, Structural Dynamics, and Materials 
Conf., April 20-23, 1998, Long Beach, CA. 


8 


American Institute of Aeronautics and Astronautics 



Figures 


Spray uncertainty onto inside of master structure . . . 




Figure 1: A FE model is made fuzzy by attaching appropriate boundary conditions to nodes. Each 

boundary condition is modeled as the impedance Z of an attached 1-DOF oscillator. Note the the total 
number of DOFs in the model is unchanged when the fuzzy substructure is added. 



x(t) 


Figure 2: The 1-DOF model used for each fuzzy DOF of the fuzzy substructure. The moving support, or 

master structure, moves with displacement x(t) and velocity x(t) = u(t) = ve tuJt . The resulting force due to 
the fuzzy 1-DOF attachment is F(t) = Fe tUit . 
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-40 


master structure (dashed), with fuzzy (solid) 



Figure 3: Mean squared acceleration at a point on Soize’s simply supported beam versus frequency. The 

dashed line shows the result for the master structure alone, while the solid line shows the result including 
the effects of the fuzzy substructure. The horizontal axis is the frequency in hertz and the vertical axis is the 
expected mean squared acceleration in dB relative to 1 m/s 2 . A 1 Hz frequency resolution is used. These 
results compare favorably to Fig. 7(a) of Soize’s work. 5 
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